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Abstract — Consider a network where all nodes are distributed 
on a unit square following a Poisson distribution with known 
density p and a pair of nodes separated by an Euclidean 
distance x are directly connected with probability g where 

g : [0, oo) — > [0, 1] satisfies three conditions: rotational in- 
variance, non-inc reasing monotonicity and integral boundedness, 

r p ~ \J '° e cJ' b ' C ~ fm 9 (W X W) and 6 is a constant, 
independent of the event that another pair of nodes are directly 
connected. In this paper, we analyze the asymptotic distribution 
of the number of isolated nodes in the above network using the 
Chen-Stein technique and the impact of the boundary effect on 
the number of isolated nodes as p — > oo. On that basis we derive 
a necessary condition for the above network to be asymptotically 
almost surely connected. These results form an important link 
in expanding recent results on the connectivity of the random 
geometric graphs from the commonly used unit disk model to 
the more generic and more practical random connection model. 



Index Terms- 
model 



-Isolated nodes, connectivity, random connection 



I. Introduction 



Connectivity is one of the most fundamental properties of 
wireless multi-hop networks 0~)-||5). A network is said to 
be connected if there is a path between any pair of nodes. 
In this paper we consider the necessary condition for an 
asymptotically almost surely (a.a.s.) connected network in 5ft 2 . 
Specifically, we investigate a network where all nodes are 
distributed on a unit square [—5)5) following a Poisson 
distribution with known density p and a pair of nodes sep- 
arated by an Euclidean distance x are directly connected with 
probability g independent of the event that another pair 

of nodes are directly connected. Here g : [0, 00) — > [0, 1] 
satisfies the properties of rotational invariance, non-increasing 
monotonicity and integral boundedness [6|, [7 Chapter 6] [if 



g{x) <g (y) 

< fgp g (\\x\\) dx < 00 



whenever x > y 



(1) 



where r p - 
constant and 



log p+b 



- Up — , < C = J m g(\\x\\)dx < 00, b is a 
denotes the Euclidean norm. 



This research is funded by ARC Discovery project: DP0877562. 

'Throughout this paper, we use the non-bold symbol, e.g. x, to denote a 
scalar and the bold symbol, e.g. x, to denote a vector. 

2 We refer readers to (6), [JJ Chapter 6] for detailed discussions on the 
random connection model. 



It is shown later in Section III-BI that the conditions in 
(Q3 imply g (x) — o x \\j where the symbol o x is defined 
shortly later. In this paper we further require g to satisfy a 
slightly more restrictive condition that g (x) = o x ( x2 r^p 
and the implications of such more restrictive condition become 
clear in the analysis of Section III-BI particularly in Remark 
[TJ The condition g(x) = o x ( x2 - ^ is only slightly more 
restrictive than the condition g (x) = o x (^j-J in that for an 
arbitrarily small positive constant e, -h? = o r ( « , 1 a ) . 

j r x A ~r£ ^ I X A log^ X J 

The reason for choosing this particular form of r p is that the 
analysis becomes nontrivial when b is a constant. Other forms 
of r p can be accommodated by allowing b —> 00 or b —> —00, 
e.g. b becomes a function of p, as p — > 00. We discuss these 
situations separately in Section [IV] 

Denote the above network by Q (X p ,g p ). It is obvious that 
under a unit disk model where g(x) — 1 for x < 1 and g(x) — 
for x > 1, r p corresponds to the transmission range for 
connectivity JTJ. Thus the above model easily incorporates the 
unit disk model as a special case. A similar conclusion can 
also be drawn for the log-normal connection model. 

The following notations and definitions are used: 

. f(z) = o z (h (z)) iff (if and only if) lim^oo ^ = 0; 
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>00 l,f~\ -L? 



• / (z) ~ z h (z) iff lim^oo 

• An event £ z depending on z is said to occur a.a.s. if its 
probability tends to one as z — > 00. 

The above definition applies whether the argument z is con- 
tinuous or discrete, e.g. assuming integer values. 

The contributions of this paper are: firstly using the Chen- 
Stein technique 0, J5), we show that the distribution of 
the number of isolated nodes in Q(X p ,g p ) asymptotically 
converges to a Poisson distribution with mean er h as p —> 00; 
secondly we show that the number of isolated nodes due to the 
boundary effect in Q (X p ,g p ) is a.a.s. zero, i.e. the boundary 
effect has asymptotically vanishing impact on the number of 
isolated nodes; finally we derive the necessary condition for 
Q (X p ,g p ) to be a.a.s. connected as p — > 00 under a generic 
connection model, which includes the widely used unit disk 
model and log-normal connection model as its two special 
examples. 

The rest of the paper is organized as follows: Section [TT] 
analyzes the distribution of the number of isolated nodes on 
a torus; Section [III] evaluates the impact of the boundary 



effect on the number of isolated nodes; Section [TV] provides 
the necessary condition for Q (X p ,g p ) to be a.a.s connected; 
Section [V] reviews related work in the area. Discussions on 
the results and future work suggestions appear Section |VT] 

II. The Distribution of The Number Of Isolated 
Nodes on A Torus 

Denote by Q T (X p ,g p ) a network with the same node 
distribution and connection model as G(X p ,g p ) except that 
nodes in Q T (X p ,g p ) are distributed on a unit torus [— \, i) . 
In this section, we analyze the distribution of the number 
of isolated nodes in Q T (X p ,g p ). With minor abuse of the 
terminology, we use A T to denote both the unit torus itself 
and the area of the unit torus, and in the latter case, A T = 1. 



A. Difference between a torus and a square 



The unit torus [- 



1 ii 



that is commonly used in random 
geometric graph theory is essentially the same as a unit 
square [—5, 5) except that the distance between two points 
on a torus is defined by their toroidal distance, instead of 
Euclidean distance. Thus a pair of nodes in Q T (X p ,g p ), 
located at X\ and x^ respectively, are directly connected 
with probability g ^ H 3il ~ 3:2 ll ^ where ||a;i — 2:2 || T denotes 
the toroidal distance between the two nodes. For a unit torus 
tj, :j) , the toroidal distance is given by iflOl p. 13]: 

{\\x 1 +z-x 2 \\:zeZ 2 } (2) 



a t = r-i 



11 T A 

\x\ — X2W = mm 



The toroidal distance between points on a torus of any other 
size can be computed analogously. Such treatment allows 
nodes located near the boundary to have the same number of 
connections probabilistically as a node located near the center. 
Therefore it allows the removal of the boundary effect that is 
present in a square. The consideration of a torus implies that 
there is no need to consider special cases occurring near the 
boundary of the region and that events inside the region do 
not depend on the particular location inside the region. This 
often simplifies the analysis however. From now on, we use the 
same symbol, A, to denote a torus and a square. Whenever the 
difference between a torus and a square affects the parameter 
being discussed, we use superscript T (respectively s ) to mark 
the parameter in a torus (respectively a square). 

We note the following relation between toroidal distance and 
Euclidean distance on a square area centered at the origin: 

\T 11 11 



Xl 



*al| T < 



Ixi — X2W and 



(3) 



which will be used in the later analysis. 



B. Properties of isolated nodes on a torus 

Divide the unit torus into m 2 non-overlapping squares each 
with size -\. Denote the square by Ai m . Define two sets 
of indicator random variables Jf and If with i m E T. m = 
{1, . . . m 2 }, where Jf — 1 iff there exists exactly one node in 
A im , otherwise Jf — 0; if = 1 iff there is exactly one node 
in Ai m and that node is isolated, If = otherwise. Obviously 
Jf m is independent of Jj m ,j m € T m \ {i m }. Denote the center 



of Af by Xi m and without loss of generality we assume 
that when jf = 1, the associated node in A im is at x i7 J^. 
Observe that for any fixed m, the values of Pr (if = l) and 
Pr (jf = l) do not depend on the particular index i m on 
a torus. However both the set of indices r m and a particular 
index i m depend on m. As m changes, the square associated 
with If and Jf also changes. Without causing ambiguity, 
we drop the explicit dependence on m in our notations for 
convenience. As an easy consequence of the Poisson node 
distribution, 

Pr(Jf = l 



lim 



1 



(4) 



and as m — y 00, the probability that there is more than 
one node in A, becomes vanishingly small compared to 
Pr ( Jf = l). Further, using the relationship that 



Pr (If = 1) - Pr (If 
it can be shown that 

Pr (If = 1) 
= Pr(jf = l) 



l\jf = l) Pr(jf = l) 



(5) 



n 

ier\{i} 



PriJi 



1 1 



1-Pr J = U-Om Pr J. 



Pr (Jf = 11-5 




g f ilm jt>i\\ represents 



In ©, the term Pr (jj - 

the probability of the event that there is a node in Aj and 
that node is not directly connected to the node in Ai, the 
term (l — Pr (jj = l) — o m (Pr (jj = l))) represents the 
probability of the event that there is no node in Aj and the 
last term accounts for the situation that there is more than one 
node in Aj. It then follows that 

lim Prilf = 1\J, = l) 



= lim 



n 

jer\{i} L 



- I A P3 



1 - Pr 



dx 



W - 1) 



IaP9 



<1'X 



(7) 
(8) 

where © results from due to nodes being distributed on 
a torus. Further, using (0), (0 and (0, it is evident that 



JaP9 



tlx 



(9) 



Define W% = J2T=i 7 f and wT = lim »^oo W%, where 
W T has the meaning of the total number of isolated nodes in 
A. It then follows that 



E (W T ) = lim E (Wf) = pe 



SaPS 



dx 



(10) 



In this paper we are mainly concerned with the case that m — > 00, i.e. 
the size of the squer is varnishingly small. Therefore the actual position of 
the node in the square is not important. 



It can be shown that 



lim pe 




p— foo 




= lim pe 




p— s-oo 








= lim pe 




p— 5-00 






log P + fc 


= e~ b lim 


c 

e 



p— >oo 

where Z? (0, a;) denotes a disk centered at the origin and with 
a radius x, e is a small positive constant, and the last step 
results because 

r°? e 2irxg (x) dx 
lim — ^ 

p— >oo 



lim 

p— >oo 



logp+6 



p(log p+6)^ 

lim 7re (log p + b) 



p— >oo 



2 -2e 
r p Op 



„-2e i„„2 / -2e 



log (r-p 



(12) 



= 



where UHi^cepital's rule is used in reaching ( fTSl i and 
in the third step g (x) 



"^i^i^l is used - As 



a consequence of (O, ( fTQb , dTTb and that 



,-6 



linip^oo pe 
limp_>. 



•feO* < lim 
•/ JJ (o^-.) W (¥)« fa 



lim £ (VF r ) 

p— >-00 V y 



< 



(13) 



The above analysis is summarized Lemma Q] 

Lemma 1: The expected number of isolated nodes in 



S T (X P ,g P ) is pe 



<l'X 



As p — > oo, the expected 



number of isolated nodes in Q (X p ,g p ) converges to 



Remark 1: Using £[), it can 

C = J u2 g (\\x\\) dx > lim 2 ^ f 



be shown that 

J n 27rxg (2) cfa; = 
< -. It can 



9( z ) _ g 

9(2) 



(14) 



lirriz^oo 7rz 2 g (z). Therefore lim^ 

then be shown that the only possibility is lim 

and that the other possibilities where lim z _ i . 00 ^ can 
be ruled out by contradiction with ((TJ. Thus 

g(x) = o x (l/x 2 ) 

Further the condition g (x) = o x y x2 ^5 - ^ is only re- 
quired for pr 2 p §sft2\ D (Q r -e\ g (\\ x \\) dx to asymptotically con- 
verge to 0, where the term pr 2 p J%p\ D t Q r -<^ 9 {\\ x \\) dx is 
associated with (the removal of) connections between a 
node at and other nodes outside D(0,r~ e ). Evalua- 
tion of pr 2 J^2\ D ^ g (\\x\\) dx for an area larger than 

D (0, r p £ ) (but not greater than A p ) does not remove the need 
for the condition. Thus the more restrictive requirement on g 
that g (x) = o x ^ x2 ^ 2 - ) is attributable to the truncation 



effect that arises when considering connectivity in a (asymp- 
totically infinite) finite region instead of an infinite area. 
Now consider the event if I J — l,i 7^ j conditioned on the 
event that JfJf = 1, meaning that both nodes having been 
placed inside A; and Aj respectively are isolated. Following 
the same steps leading to ®, it can be shown that 

lim Pr (if if = 1| Jf Jf = 1) 



x exp 




where the term (l — g ( " Xt r ffijl1 j j is due to the consid- 
eration that the two nodes located inside Ai and Aj cannot 
be directly connected in order for both nodes to be isolated. 
Observe also that: 

Pr {U Ij = 1) 
= Pr (jfjf = 1) Pr (If If = l\jfjf = 1) (16) 

Now using ©, (O, (H~5T > and ( fTBI l. it can be established that 

Pr (if = l\lT = 1) 




,. ( „ (s (fc^ H t^ H fc^)) 



(17) 



lim 



Pr (IT If = 1) 



Pr (If = 1) Pr (If = 1) 



= (i _ g( m - si! ))e /x Pg (^^) g (^^)ti8) 

Using ©, ((TJ, (JT3J, ( fT6] > and the above equation, it can also 
be obtained that 



Pr 1 7 = 1,1 T = 



lim 

m->oo p r in = 1) Pr (2? 



0) 



lim 



Pr (If = l)-Pr (If If = 1 
Pr (If = 1) Pr (If = 0) 

-jir 



,. f, p ~SaP9 

lim 1 ^e 

m->oo \ m 



x e 




(19) 



C. The asymptotic distribution of the number of isolated nodes 

On the basis of the discussion in the last subsection, in 
this subsection we consider the distribution of the number 
of isolated nodes in Q T (X p ,g p ) as p — > oo . Our analysis 
relies on the use of the Chen-Stein bound JS], 0. The Chen- 
Stein bound is named after the work of Stein [ 1 1 1 and Chen 
H21, lfl3ll . It is well known that the number of occurrences of 
independently distributed rare events in a region can often 
be accurately approximated by a Poisson distribution (9). 
In ifTTI . Stein developed a novel method for showing the 
convergence in distribution to the normal of the sum of a 
number of dependent random variables. In lfl2l . Ifl3l Chen 
applied Stein's ideas in the Poisson setting and derived an 
upper bound on the total variation distance, a concept defined 
in the theorem statement below, between the distribution of 
the sum of a number of dependent random indicator variables 
and the associated Poisson distribution. The following theorem 
gives a formal statement of the Chen-Stein bound. 

Theorem 1: J9] Theorem l.A ] For a set of indicator 
random variables Ii, i e T, define W = Xier Pi — E {h) 
and A = E(W). For any choice of the index set T St i C I\ 
r s ,i n {i} = -foce}, 

d TV {C{W),Po{\)) 



< 



E 




+ H E 7 * E Ij)min(l,j 

ier \ jer s ,t j x 

+ e E \ E ^ i&> j e r «.*) y-p*\ mm ( X < ^) 

where C (W) denotes the distribution of W, Po (A) denotes a 
Poisson distribution with mean A, T W: i — Y\ {T s ^ U {i}} and 
dxv denotes the total variation distance. The total variation 
distance between two probability distributions a and f3 on Z + 
is defined by 

d TV { a , fi) = sup {\a {A) -(3{A) \ :AcZ+} 

For convenience, we separate the bound in Theorem [TJ into 
three terms b\min (l, jj, bivnin (l, j) and b^min ^1 



where b\ = 
The set of indices T, 



(/*E : 



pi 



i 



Pi- 5 i t^jer, 
and 63 



/ Pi]' 

is often chosen to contain all those 



j, other than i, for which Ij is "strongly" dependent on Ii 
and the set often contains all other indices apart from 
i for which Ij,j £ Tw,i are at most "weakly" dependent on 
h El. In many applications, by a suitable choice of T s ,i the 
&3 term can be easily made to be and the evaluation of 
the b\ and &2 terms involve the computation of the first two 
moments of W only, which can often be achieved relatively 
easily. An example is a random geometric network under the 
unit disk model. If T Sy i is chosen to be a neighborhood of i 
containing indices of all nodes whose distance to Xi is less 



than or equal to twice the transmission range, the 63 term is 
easily shown to be 0. It can then be readily shown that the b\ 
and 62 terms approach as the neighbourhood size of a node 
becomes vanishingly small compared to the overall network 
size asp->oo [14|. However this is certainly not the case for 
the generic random connection model where the dependence 
structure is global. 

Using the Chen-Stein bound, the following theorem, which 
summarizes a major result of this paper can be obtained: 

Theorem 2: The number of isolated nodes in Q T {X p ,g p ) 
converges to a Poisson distribution with mean e~ b as p — > 00. 
Proof: Proof is given in the Appendix. ■ 

III. The Impact of the Boundary Effects on the 
Number of Isolated Nodes 

On the basis of the analysis in the last section, we now 
consider the impact of the boundary effect on the number of 
isolated nodes in Q (X p ,g p ). Following the same procedure 
that results in ©, it can be shown that Pr (if = l) ~ m 

_p_ e /a r-p \\) dx w here the parameters in this section 
is defined analogously as those in the last section. Note that 
due to the consideration of a square, a relationship such as 
JaP9 ( 7~ ) ^ x is no longer valid. 

Ja«>(|| 2 



/a«( 



It follows that 



dx 



E (W s ) = Km E (W£) = 



pe 



) dx 



dy 



lim E (W s 

p—>oo 



= lim 

p—>-oo 



S An pr 2 p g{\\x-y\\)dx 



lim pe 

p— J- 00 



-Cprt 



dy 
(20) 



2b 2H - In 



where A p is a square of size 3 and A p = 
arriving at < f20b some discussions involving dividing A p into 
three non-overlapping regions: four square areas of size r~ e x 
x p , denoted by ZA p ; four rectangular 



t.— e 

' p 



areas of size r p e x [r~ — 2r p £ ) adjacent to the four sides 
of A p , denoted by £A p ; and the rest central area, are omitted 
due to space limitation, where e is a small positive constant 
and e < j. 

Comparing (fTJt and (|20| >. it is noted that the expected num- 
bers of isolated nodes on a torus and on a square respectively 
asymptotically converge to the same nonzero finite constant 
e~ b as p —> 00. Now we use the coupling technique I0 to 
construct the connection between W s and W T . Consider an 
instance of Q T (X p ,g p ). The number of isolated nodes in that 
network is W T , which depends on p. Remove each connection 

a(f-) 

of the above network with probability 1 > , independent 

3{Tp-) 

of the event that another connection is removed, where x 
is the Euclidean distance between the two endpoints of the 
connection and x T is the corresponding toroidal distance. Due 

g(f-) 

to (O and the non-increasing property of q, < 1 } §A. < 

1. Further note that only connections between nodes near 
the boundary with x T < x will be affected. Denote the 
number of newly appeared isolated nodes by W E . W E has 



the meaning of being the number of isolated nodes due to the 
boundary effect. It is straightforward to show that W is a 
non-negative random integer, depending on p. Further, such 
a connection removal process results in a random network 
with nodes Poissonly distributed with density p where a pair 
of nodes separated by an Euclidean distance x are directly 
connected with probability g (^-)> i- e - a random network on 
a square with the boundary effect included. The following 
equation results as a consequence of the above discussion: 
W s = W E + W T . Using (QI), (|20]l and the above equation, 
it can be shown that lim^oo E (W E ) = 0. Due to the non- 
negativeness of W E : lirrip^oo Pr [W E = 0) = 1. The above 
discussion is summarized in the following lemma, which forms 
the second major contribution of this paper. 

Theorem 3: The number of isolated nodes in Q (X p , g p ) due 
to the boundary effect is a.a.s. as p — > oo. 

IV. The Necessary Condition for Asymptotically 
Connected Networks 

We are now ready to present the necessary condition for 
Q (X p ,g p ) to be a.a.s. connected as p — > oo. The following 
theorem can be obtained using Theorems [2] and [3] 

Theorem 4: The number of isolated nodes in Q(X p ,g p ) 
converges to a Poisson distribution with mean e~ b as p — > oo. 
Corollary [T] follows immediately from Theorem [4] 

Corollary 1: As p — > oo, the probability that there is no 
isolated node in Q (X p ,g p ) converges to e~ e . 
With a slight modification of the proof of Theorem [2] it can be 
shown that Theorems [2] and [4] and Corollary[T]can be extended 
to the situation when b is a function of p and lim p _j. 00 b = B, 
where B is a constant. Now we further relax the condition in 
Theorem [2] on b and consider the situation when b — > — oo 
or b — > oo as p —> oo. When ii -> oo, the number of 
connections in Q(X p ,g p ) increases. Unsurprisingly, isolated 
nodes disappear. In fact, using the coupling technique, Lemma 
[1] Theorem [3] and Markov's inequality, it can be shown that 
if b -> oo as p -> oo, lim^oo Pr {W s = 0) = 1. 

Now we consider the situation when b — > — oo as p ^ oo. 
For an arbitrary network, a particular property is increasing 
if the property is preserved when more connections (edges) 
are added into the network. A property is decreasing if 
its complement is increasing, or equivalently a decreasing 
property is preserved when connections (edges) are removed 
from the network. It follows that the property that the network 
Q (X p ,g p ) has at least one isolated node, denoted by A, is 
a decreasing property. The complement of A, denoted by 
A c , viz. the property that the network Q(X p ,g p ) has no 
isolated node, is an increasing property. In fact the network 
Qi (X p ,g p ) where b — B\ can be obtained from the network 
Qi (X p ,g p ) where b = B-i and B^ < B\ by removing each 
connection in Gi{X p ,g p ) independently with a probability 



9 



I log p+B 2 

V — e> — 



/9 



with x being the distance 



between two endpoints of the connection. The above obser- 
vations, together with Corollary [T] lead to the conclusion that 

if 6 — >• — oo as p —¥ oo, 



lim Pr(A) 

p— >oo 



lim 1 

p— >oo 



Pr (A c 



1 



The above discussions are summarized in the following 
theorem and corollary, which form the third major contribution 
of this paper: 

Theorem 5: In the network Q (X p ,g p ), if b — ► oo as p — > 
oo, a.a.s. there is no isolated node in the network; if b — > — oo 
as p — > oo, a.a.s. the network has at least one isolated nodes. 

Corollary 2: b — > oo is a necessary condition for the 
network Q (X p ,g p ) to be a.a.s. connected as p — > oo. 

V. Related Work 

Extensive research has been done on connectivity problems 
using the well-known random geometric graph and the unit 
disk model, which is usually obtained by randomly and uni- 
formly distributing n vertices in a given area and connecting 
any two vertices iff their distance is smaller than or equal to 
a given threshold r(n) [10], [15|. Significant outcomes have 
been achieved for both asymptotically infinite n JT], J2), [ 10 1, 
ED-lEO) and finite n 01. HI. EH. Specifically, it was shown 
that under the un it disk mo del and in 5ft 2 , the above network 

with r (n) = J log n ^f ^- is a.a.s. connected as n oo 
iff c(n) — !• oo . In [17|, Ravelomanana investigated the 
critical transmission range for connectivity in 3-dimensional 
wireless sensor networks and derived similar results as the 2- 
dimensional results in |1|. Note that most of the results for 
finite n are empirical results. 

In 0, Il22l - ll26l the necessary condition for the above 
network to be asymptotically connected is investigated under 
the more realistic log-normal connection model. Under the log- 
normal connection model, two nodes are directly connected if 
the received power at one node from the other node, whose 
attenuation follows the log-normal model, is greater than a 
given threshold. These analysis however all relies on the 
assumption that the event that a node is isolated and the event 
that another node is isolated are independent. Realistically 
however, one may expect the above two events to be correlated 
whenever there is a non-zero probability that a third node may 
exist which may have direct connections to both nodes. In 
the unit disk model, this may happen when the transmission 
range of the two nodes overlaps. In the log-normal model, 
any node may have a non-zero probability of having direct 
connections to both nodes. This observation and the lack of 
rigorous analysis on the node isolation events to support the 
independence assumption raised a question mark over the 
validity of the results of 0, Il22l-ll26l. 

The results in this paper complement the above studies in 
two ways. They provide the asymptotic distribution of the 
number of isolated nodes in the network, which is valid not 
only for the unit disk model and the log-normal connection 
model but also for the more generic random connection model. 
Second they do not depend on the independence assumption 
concerning isolated nodes just mentioned. In fact, it is an 
unjustifiable assumption. They do however rely on the inde- 
pendence of connections of different node pairs, referred to in 
the discussion of the random connection model in Section U 

Some work exists on the analysis of the asymptotic distri- 
bution of the number of isolated nodes 0, HSl, [3, |27l 
under the assumption of a unit disk model. In ll27ll . Yi et 



al. considered a total of n nodes distributed independently 
and uniformly in a unit-area disk. Using some complicated 
geometric analysis, they showed that if all nodes have a 
maximum transmission range r(n) = y (log n + £) /7m for 
some constant £, the total number of isolated nodes is asymp- 
totically Poissonly distributed with mean e~^. In (6), [14|, 
Franceschetti et al. derived essentially the same result using 
the Chen-Stein technique. A similar result can also be found 
in iflOl in a continuum percolation setting. There is a major 
challenge in analyzing the distribution of the number of 
isolated nodes under the random connection model; under the 
unit disk model, the dependence structure is "local", i.e. the 
event that a node is isolated and the event that another node is 
isolated are dependent iff the distance between the two nodes 
is smaller than twice the transmission range, whereas under the 
random connection model, the dependence structure becomes 
"global", i.e. the above two events are dependent even if the 
two nodes are far away. 

VI. Conclusions and Further Work 

In this paper, we analyzed the asymptotic distribution of 
the number of isolated nodes in Q (X p ,g p ) using the Chen- 
Stein technique, the impact of the boundary effect on the 
number of isolated nodes and on that basis the necessary 
condition for Q (X p ,g p ) to be a.a.s. connected as p — > oo. 
Considering one instance of such a network and expanding the 
distances between all pairs of nodes by a factor of l/r p while 
maintaining their connections, there results a random network 
with nodes Poissonly distributed on a square of size 1/r 2 with 
density pr 2 where a pair of nodes separated by an Euclidean 
distance x are directly connected with probability g (x). Using 
the scaling technique 0, it can be readily shown that our 
result applies to this random network. By proper scaling or 
slight modifications of the proof of Theorem [2] our result can 
be extended to networks of other sizes. 

It can be easily shown that as p — > oo, the average node 
degree in Q (X p ,g p ) converges to log p+b. That is, the average 
node degree under the random connection model increases 
at the same rate as the average node degree required for a 
connected network under the unit disk model as p — >• oo 
|16|. Further if b — >• oo as p — > oo, a.a.s. there is no 
isolated node in the network. This result coincides with the 
result in [1| on the critical transmission range required for an 
a.a.s. connected network. Another implication of our result 
is that different channel models appear to play little role 
in determining the asymptotic distribution of isolated nodes 
(hence the connectivity) so long as they achieve the same 
average node degree under the same node density. 

This paper focuses on a necessary condition for Q (X p ,g p ) 
to be a.a.s. connected. We expect that as p —> oo, the 
necessary condition also becomes sufficient, i.e. the network 
becomes connected when the last isolated node disappears. 
It is part of our future work to investigate the sufficient 
condition for asymptotically connected networks under the 
random connection model and validate the above conjecture. 

This paper focuses on the asymptotic distribution of the 
number of isolated nodes, i.e. the number of nodes with a 



node degree k = 0. We conjecture that for a generic k, the 
asymptotic distribution of the number of nodes with degree 
k may also converge to a Poisson distribution. Thus, it is 
another direction of our future work to examine the asymptotic 
distribution of the number of nodes with degree k, where 
k > 0. 

Appendix: Proof of Theorem [2] 

In this appendix, we give a proof of Theorem |2] using the 
Chen-Stein bound in Theorem Q] The key idea involved using 
Theorem[T]to prove Theorem[2]is constructing a neighborhood 
of a node, i.e. T Si i in Theorem [T] such that a) the size of the 
neighborhood becomes vanishingly small compared with A as 
p —> oo. This is required for the b\ and 62 terms to approach 
as p — > 00; b) a.a.s. the neighborhood contains all nodes that 
may have a direct connection with the node. This is required 
for the 63 term to approach as p — > 00. Such a neighborhood 
is defined in the next paragraph. 

First note that parameter W in Theorem Q] has the same 
meaning of defined in Section [TTJ Therefore the parameter 
A in Theorem [T] which depends on both p and to, satisfies 
limp-j.oo lim„ woo A = e~ b . Further Pi = E (if) and E (if) 
has been given in ([9JI. Unless otherwise specified, these pa- 
rameters, e.g. Xi, m, if ', r and r p , have the same meaning as 
those defined in Section|II] Denote by D (xi, r) a disk centered 
at Xi and with a radius r. Further define the neighbourhood 
of an index i G T as r s ^ = {j : Xj € D (a;,, 2r 1 p ~' L ) } 
and define the non-neighbourhood of the index i as T w j = 
{j : Xj £ D (a;,, 2r p ~ e ) } where e is a constant and e £ 
(0, It can be shown that 

= mHirr 2 - 2 * + o m (to 2 4^~ 2£ ) (21) 
From (O, pi = E (if) and (O, it follows that 

o — f a pq( — — — — — 1 dx 

lim m 2 pi = pe A V 1 " ) (22) 

m— >-oo 

lim lim rr?pi = e~ b (23) 

p— >oo m— >oo 

Next we shall evaluate the b\, 62 and 63 terms separately. 

A. An Evaluation of the b\ Term 
It can be shown that 

= lim m 2 Pl E V if] 
yer,, i u{i} / 

= lim (m 2 pi) 2 4irr 2 ~ 2e 

= 4*^* V J j {-^) (24) 

where in the second step, (1211 1 is used and in the final step 
(01, (l22l and the value of r p are used. It follows that 

lim RHS of m = ^ 2b lim ( log P + b ) ' = 



where ( fTOb and ( fT3l > are used in the above equation, and RHS 
is short for the right hand side. This leads to the conclusion 
that linip^oo lim m _>. 00 b\ = 0. 

B. An Evaluation of the bi Term 
For the b^ term, we observe that 
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where A p ~ 



^-,2r) , in the first step, ©, (Q3]l and 
( fT6] l are used and the final step involves some translation and 
scaling operations. Let A = los ^ +fc , it can be further shown 
that as p —> 00, 



e 2b lim RHS of (|25]l 
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logp 
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logp(logp + 6- 1) 
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lo SPl3j2h(x,y) dx 
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(26) 



where it is a unit vector pointing to the +x direction and 



x-u 



, in the first step 



log p+b 
Cp 



and 1 — 3 (j|y|| T ^ < 1 are used, 



h(x,y) = 

(Hi, (E), 

and in the last step, L'Hi^cepital's rule, where Cp is used as 
the denominator and the other terms are used as the numerator, 
(O and the following formulas are used: 

fh(x) 



d f (x> 

Txj f{x > y)dy 

f h(x) df(x,y) 
Jo 



dx 



dy + f(x,h(x)) 



dh (x) 



dx 
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, log p + b — 1 



In the following we show that all three terms inside the 



limbec sign and separated by + sign in 
p^co. First it can be shown that 



approach as 



/ 9 (\\x\\) g (\\x - u2r e \\) dx 
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(27) 



where in the second step the observation that the distance 
between any point in Z?(0,r~ e ) and u2r~ t is larger than 
or equal to r~ e , the observation that the distance between any 
point in 3i 2 \D (0, r^ 6 ) and the origin is larger than or equal 
to r~ e and the non-increasing property of g are used, (fl4] > is 
used in the last step. This readily leads to the result that the 
first term in d26b satisfies: 
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where L'Hi^oepital's rule, where Cp is used as the denom- 
inator and the other terms are used as the numerator, and 
r p = ^/ l0 g P p +b are used in the first step of the above equation, 
in the second step d27b is used, which readily leads to the 
conclusion that 



lim e 

p— >oo 



! ^J K 2 9(l|x|r) 3 (||x- U 2r p -|| T )dx = j 



(28) 



The final steps are complete by putting the value of r p into 
the equation and noting that J^ 2 h(x,y) dx < C for y ^ 0, 
which is a consequence of the following derivations: 



R 2 



(Hf) 
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dec - C 



- 1 dec < 



and the only possibility for J„j 2 /i (a;, y) dx — C = to occur 
is when g corresponds to a unit disk model and y = 0. 
For the second term in d26l l, it can be shown that 
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p— foo 
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where in the first step 
For the third term in 



is used, 
it is observed that 



= 



lim 



logp 
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h (x, y) dx 
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logp 
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where JLj 2 h (x, y) dx < C for y ^ is used in the first step. 
Eventually we get lim p _ i . 00 lim m _ i . 00 62 = 0. 

C. An Evaluation of the 63 Term 

Denote by I\ a random set of indices containing all indices 
j where j € T w ,i and Ij = 1, i.e. the node in question is 
also isolated, and denote by 7^ an instance of T^. Define n = 
Following a similar procedure that leads to ( TT~8T > and dT9l 
and using the result that j A pg ( ^ x ~™ % ^ j g ( ^ x ™ 3 - j dec = 

o p (1) and 5 ( l|x \. p Xjl|T ) = o p (1) for \\ Xi - x 3 f > 2r p ^ 
(see d28Ti). it can be shown that 
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Note that a;, and Xj,j e j is separated by a distance not 
smaller than 2r~ e . A lower bound on the value inside the 



expectation operator in d29j l is given by 



B, 



(30) 

An upper bound on the value inside the expectation operator 
in ( 129b is given by 



IaP9 



Bu,i = e 
It can be shown that 



n, 
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(32) 



Let us consider E\E{I{ \ € ^ w ,i)} ~ Pi\ now - From 
d30T >. (SB and ([32), it is clear that 

lim lim S2E\E{jT \(lf,j£r Wti )}- Pi \ 

p— >oo m—t-oo * — ' 1 1 

0, max < lim lim m 2 pi — pE (Blj) , 
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p— >-oo m— »-oo ' 



(33) 



In the following we will show that both terms 

lim TO _ J . 00 m 2 p, i - pE (B Lji ) and liirim^co pi? (B Uti ) - m 2 pi 
in ( F3~3l ) approach as p —> 00. First it can be shown that 



lim pE(B LA ) 



> lim pE (1 - ng (2r~ e )) e 



~ /a P9 



= lim p(l-S(n).g(2r; e ))e 



dec 



fix 

(34) 



where lim m _ i . 00 i? (71) is the expected number of isolated 
nodes in A\D (xi, 2rp _e ). In the first step of the above 
equation, the inequality (1 — x) > 1 — nx for < x < 1 
and n > is used. When p — > 00, rj -6 — > and r p £ — > 00 
therefore limp-j.oo lim m _ i . 00 E (n) = linip^oo E (W) = e~ b is 
a bounded value and lim^oo lim m _ s . 00 g (2r~ £ ) — > 0, which 
is an immediate outcome of ( TT~4T > . It then follows that 

lim lim — — ^ L,% ^ 

p— >oo m— >oo m pi 

> lim lim (l - E (n) g (2r~ e )) = 1 

Together with (|23]l and([32]), it follows that 

lim lim m 2 Pl - pE (B Li ) = (35) 

p— ^00 m— foo ' 

Now let us consider the second term limm^oo pE (Bjj,i) — 
m 2 pi, it can be observed that 
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(36) 



where in the second step, the non-increasing property of g, 
and the fact that Xj is located in A\D (xi,2r p ~ e ) and x is 
located in D (xi,r p ~~ e ^), therefore ||a; — Xj\\ T > r p ~ e is used. 
It can be further demonstrated, using similar steps that result 
in ([Toll and £TJ), that the term J^, r i_^ pg M |a! ~ a! ' 1 ^ da; in 
have the following property: 



X Xj t 



dx 



< Cpr 2 p = \ogp + b 











X 


1 







dx 



(37) 



For the other term ng (r~ e ) in ( |36l l, choosing a positive 
constant 6 < 2e and using Markov's inequality, it can be 
shown that 

Pr (n > r~ 5 ) ^ r s p E (n) 
lim lim Pr (ng (r~ e ) n (e, p) > r p s g (r p e ) n (e, p)) 

p^-oo m— too ' \ f / / 



< lim Jim r s p E (n) 



p—t-oo m— >-oo 



where liirip^co r p s g (r e ) n (e, p) — due to (1T4) . (|37) 

° - for any positive con- 

5, 



and 5 < 2e, linip_ s . 00 r p 



stant B, and lirrip^co lim„ l _ i . 00 r° p E (n) 
linip^oo limTO^oo E (n) = lim p _>. 00 £7 (W) 
bounded value and that lim p _ i . 00 r p = 0. Therefore 



due to that 

e~ b is a 



lim lim Pr (ng (r e ) r; (e, p) = 0) = 1 (38) 

As a result of (06), (J38>, (02), flO) and (E): 

lim lim pP (Bui) 

p— >oo m— ^oo ' 



< lim lim pP e 



dec 



lim pe' 

p— >oo 



-Cprt 



Using the above equation, ( 123) andd32), it can be shown that 

lim lim pE (B Utl ) - m 2 Pl = (39) 

p— foo m— feso ' 

As a result of (133) . d35) and (|39) , limp-j.oo limm^oo 6 3 = 0. 
A combination of the analysis in subsections IVI-AI IVI-BI 
and IVI-Cl completes this proof. 
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